Angiogenesis regulators S100A4, SPARC and SPP1 correlate with macrophage infiltration and are prognostic biomarkers in colon and rectal cancers

Introduction Increasing evidence suggests that it is necessary to find effective and robust clinically validated prognostic biomarkers that can identify “high-risk” colorectal cancer (CRC) patients. Currently, available prognostic factors largely include clinical-pathological parameters and focus on the cancer stage at the time of diagnosis. Among cells of tumor microenvironment (TME) only Immunoscore classifier based on T lymphocytes showed high predictive value. Methods In the present study, we performed the complex analysis of mRNA and protein expression of crucial regulators of tumor angiogenesis and tumor progression, expressed by tumor-associated macrophages (TAMs): S100A4, SPP1 and SPARC. Colon and rectal cancer patients were investigated independently and in a combined cohort (CRC). For mRNA expression, we analyzed RNA sequencing data obtained from TCGA (N=417) and GEO (N=92) cohorts of colorectal cancer patients. For protein expression, we performed IHC digital quantification of tumor tissues obtained from 197 patients with CRC treated in the Department of abdominal oncology in Clinics of Tomsk NRMC. Results High S100A4 mRNA expression accurately predicted poor survival for patients with CRC independently of cancer type. SPARC mRNA level was independent prognostic factors for survival in colon but not in rectal cancer. SPP1 mRNA level had significant predictive value for survival in both rectal and colon cancers. Analysis of human CRC tissues revealed that S100A4, SPP1 and SPARC are expressed by stromal compartments, in particular by TAMs, and have a strong correlation with macrophage infiltration. Finally, our results indicate that chemotherapy-based treatment can change the predictive direction of S100A4 for rectal cancer patients. We found that S100A4 stromal levels were higher in patients with better response to neoadjuvant chemotherapy/chemoradiotherapy, and S100A4 mRNA levels predicted better DFS among non-responders. Discussion These findings can help improve the prognosis of patients with CRC based on S100A4, SPP1 and SPARC expression levels.

Introduction: Increasing evidence suggests that it is necessary to find effective and robust clinically validated prognostic biomarkers that can identify "high-risk" colorectal cancer (CRC) patients. Currently, available prognostic factors largely include clinical-pathological parameters and focus on the cancer stage at the time of diagnosis. Among cells of tumor microenvironment (TME) only Immunoscore classifier based on T lymphocytes showed high predictive value.
Methods: In the present study, we performed the complex analysis of mRNA and protein expression of crucial regulators of tumor angiogenesis and tumor progression, expressed by tumor-associated macrophages (TAMs): S100A4, SPP1 and SPARC. Colon and rectal cancer patients were investigated independently and in a combined cohort (CRC). For mRNA expression, we analyzed RNA sequencing data obtained from TCGA (N=417) and GEO (N=92) cohorts of colorectal cancer patients. For protein expression, we performed IHC digital quantification of tumor tissues obtained from 197 patients with CRC treated in the Department of abdominal oncology in Clinics of Tomsk NRMC.
Results: High S100A4 mRNA expression accurately predicted poor survival for patients with CRC independently of cancer type. SPARC mRNA level was independent prognostic factors for survival in colon but not in rectal cancer. SPP1 mRNA level had significant predictive value for survival in both rectal and colon cancers. Analysis of human CRC tissues revealed that S100A4, SPP1 and SPARC are expressed by stromal compartments, in particular by TAMs, and have a strong correlation with macrophage infiltration. Finally, our results indicate that chemotherapy-based treatment can change the predictive direction of S100A4 for rectal cancer patients. We found that S100A4 stromal levels were higher in patients 1 Introduction Colorectal cancer (CRC) is the third most common malignancy and the second-leading cause of cancer death in the world due to the unmet screening programs, therapeutic strategies, and increasing incidence rates (1). Colorectal cancer is characterized by high inter-and intra-tumoral heterogeneity (2,3). Although colorectal cancer is a more general and widely used term, there is still a differentiation into two distinct localizations: colon cancer (CC) and rectal cancer (RC). There are evidences accumulated towards considering CC and RC as self-standing tumor entities due to their topography, surgical challenge, therapy, complications, and relapse patterns (4,5). CRC heterogeneity determines difficulties in choosing anticancer treatment, and also poses an obstacle in reaching therapeutic complete response (6). Although the response rate to systemic chemotherapies goes up to 50%, nearly all patients with CRC develop drug resistance, which limits the therapeutic efficacy of anticancer drugs and ultimately leads to chemotherapy failure (7,8). These difficulties pose a demand for biomarker discovery that will help in improving treatment efficiency, early detection, and could be of value as diagnostic or prognostic markers.
Tumor microenvironment (TME), which consists of stromal and immune cells, has an essential role in tumor development (9). The key component of innate immunity in the TME is tumorassociated macrophages (TAMs) (10). TAMs regulate tumor growth by supporting cancer cell survival and proliferation, angiogenesis, and metastasis, as well as the response of cancer cells to therapeutic intervention (11). Angiogenesis is a basic process that provides the tumor with crucial nutrients and oxygen (12). The main pro-angiogenic regulator in tumors is VEGF (12). Despite the growing list of FDA-approved anti-VEGF drugs, the success of anti-angiogenic therapy is limited. Failure in VEGF-targeted therapy can be explained by the switching on the alternative pro-angiogenic activators (13).
Recently collected data demonstrated that TAMs can be essential sources for plenty of novel angiogenesis-related proteins that belong to S100A class, SEMA family, chitinase-like proteins, growth factors, and proteins regulating cell-matrix interactions, etc (13). Among them, pro-angiogenic factors S100A4 and osteopontin (OPN, or SPP1) as well as anti-angiogenic factor SPARC, which have drawn our attention, since we have recently identified the deregulation of their expression in TAMs under chemotherapy exposure in vitro (unpublished data).
In the present study, for the first time we performed the complex analysis of mRNA and protein expression of S100A4, SPP1 and SPARC. We establish their prognostic significance in terms of survival rates and clinical and pathological parameters of tumor state in patients with colon and rectal cancer. Also, we identified that neoadjuvant chemotherapy/chemoradiotherapy can reverse their predictive value in more favorable way.

Dataset analysis
The Cancer Genome Atlas (TCGA) data (TCGA-COAD and TCGA-READ datasets) and NCBI GEO (GSE190826 dataset) were used to examine the expression of SPP1, S100A4 and SPARC and to perform survival analysis in colorectal cancer patients. TCGA data included information about SPP1, S100A4 and SPARC expression, that was evaluated in the following groups of patients: a) with colorectal cancer (common group) (N=417), b) with colon cancer, including transverse colon, ascending colon, descending colon, sigmoid colon, cecum, hepatic flexure, splenic flexure (n=305), c) with rectal cancer, including rectosigmoid junction and rectum (N=112), with available clinical information and records on recurrence and survival rates (in details in Supplementary Table  S1). Patients with advanced stage IV were excluded. GSE190826 dataset included 92 patients with rectal cancer treated with neoadjuvant chemoradiotherapy (NCRT); information about pretreatment levels of SPP1, S100A4 and SPARC mRNA expression was obtained. The TCGA biolinks was used for retrieving RNA-seq data from the GDC database. The raw sequencing reads were processed via the DESeq2 R package. The raw counts were depth normalized and variance stabilized via the variance stabilizing transformation (VST) for downstream survival analysis.

CIBERSORT analysis
Cell type deconvolution of colorectal cancer TCGA RNA-seq data was performed via the TIMER2.0 (14) web platform, which provides a facility for robust estimation of immune infiltration levels of userprovided tumor profiles. TIMER2.0 utilizes the immunedeconv (15), an R package which integrates six distinct cell-type deconvolution algorithms, including CIBERSORT (16). CIBERSORT uses highly robust-to-noise linear support vector regression (SVR) to deconvolve the mixture of cell types of interest. Inferred immune cell-types were used to assess cell-type association of SPP1, S100A4, and SPARC genes by the Spearman correlation.

Clinical material
The IHC study included patients with colorectal adenocarcinoma with morphologically verified diagnosis, treated in the Department of abdominal oncology, Cancer Research Institute of Tomsk National Research Medical Center (Tomsk, Russia). The study was carried out according to Declaration of Helsinki (from 1964, revised in 1975 and 1983) and was approved by the local committee of Medical Ethics of Tomsk Cancer Research Institute; all patients signed informed consent for the study. Patients were divided as we did for TCGA cohort with the exception that the number of patients in SPARC/ SPP1 group differed from S100A4 group: a) with colorectal cancer (common group for SPARC/SPP1) (N=118), b) with colon cancer (N=54), c) with rectal cancer (N=64) (Supplementary Table S2). For S100A4 group: a) with colorectal cancer (common group) (N=197), b) with colon cancer (N=89), c) with rectal cancer (N=107) (Supplementary Table S3). Patients with rectal cancer and cancer of the rectosigmoid junction received neoadjuvant chemotherapy (NAC) or chemoradiotherapy (NCRT). Five-grade Mandard Tumor Regression Grading (TRG) system was used for assessment of response in patients, where TRG1no residual cancer, TRG2residual isolated cancer cells, TRG3 -fibrosis outgrowing residual cancer, TRG4residual cancer outgrowing fibrosis, TRG5absence of regressive changes (14). All patients underwent surgical treatment. In adjuvant regime, according to indications, patients received chemotherapy under the same schemes for up to 6 months. Cases of stage IV disease were excluded.

Digital quantification
Tumor tissue slides were scanned by using the Leica Aperio AT2 histoscanning station (Leica, Germany) and ScanScope software (Aperio ScanScope XT Leica). QuPath software (free from https:// qupath.github.io) was used to analyze and quantify marker expression. Individual tumor regions were selected and analyzed using cell detection and cell intensity classification. "Cell: DAB OD mean" was used for the analysis of both membranous and cytoplasmic staining of SPP1, S100A4 and SPARC. Intensity thresholds were set to further subclassify cells as being negative, weak (1+), moderate (2+) or strongly positive (3+) for marker staining based upon mean nuclear DAB optical density. The results were obtained in two scales: percentage of positive cells among all counted cells per section (%) and H-scorea parameter that takes into account both the percentage of positive cells and the intensity of staining. These parameters were counted automatically by the program.

NGS-GeoMx Digital Spatial Profiler (DSP) analysis
NanoString GeoMx digital spatial profiling (DSP) was applied to perform spatially resolved RNA profiling analysis in colorectal cancer tissue. The Cancer Transcriptome Atlas (CTA) panel was used. The 97 areas of illumination (AOIs) across all slides in mixed stroma/tumor regions was selected. The resulted libraries were sequenced by the Illumina NextSeq 500 platform using 2 x 27 base paired reads. The raw counts were processed in the NanoString's GeoMx NGS pipeline v.2.1 where they were converted to the digital count conversion (DCC) files. The GeomxTools was used for quality control (QC) and downstream analysis of the DCC files in R (17). The adjusted p-values were calculated using the Benjamini-Hochberg correction. Differential gene expression analysis between CK+ and CD45+ regions was performed using a linear mixed model (LMM) with random slope and random intercept as recommended in the GeomxTools manual.

Statistical analysis
Statistical analysis was performed using STATISTICA 12.0 for Windows (STATISTICA, RRID : SCR_014213) and GraphPad Prism 8.4.2 (GraphPad Prism, RRID : SCR_002798). The Mann-Whitney Utest and t-test for independent groups were implemented. The prognostic values of SPP1, S100A4 and SPARC (area under curve (AUC), confidence interval (CI), sensitivity, specificity, and cut-off value) were determined using receiver operating characteristic (ROC) analysis. The survival rates were determined by the Kaplan-Meier method, and the log-rank test was used. Cox's proportional-hazard model was applied for survival analysis and the hazard ratio (HR [95%CI]) evaluation. Results were presented using GraphPad Prism 8.4.2 software. Results were considered to be significant with p<0,05. Data with marginal significance (p-value >0,05 and <0,1) were also discussed.

Results
3.1 High S100A4 mRNA expression is a robust predictor for shorter survival rates independent of cancer type S100A4 is a pro-angiogenic factor belonging to the family of calcium-binding proteins (13). In the pathogenesis of cancer, increased S100A4 expression correlates with a high incidence of metastasis and poor prognosis in cancer (18).
To reveal the prognostic value of S100A4 mRNA expression, we performed survival analysis using TCGA data. The mRNA expression was defined further as log-normalized counts. We analyzed a combined group of CRC patients and two separate cohorts of patients with colon cancer and rectal cancer.
We applied Cox regression analysis, ROC analysis, and Kaplan-Meier curves to evaluate significance of S100A4 on disease prognosis in the combined CRC group, as well as CC and RC patients. In each group, all patients were categorized into high-and low-risk groups based on S100A4 expression according to the cut-off meanings determined by ROC analysis that allowed to predict overall death, death from the disease, recurrence and progression (Supplementary Table S4). Kaplan−Meier survival curve indicated that CRC patients in S100A4 high-risk group (>cut-off) had shorter overall survival (OS) (p=0,0159), disease-specific survival (DSS) (p<0,0001), disease-free survival (DFS) (p<0,0001), and progression-free survival (PFS) (p<0,0001) than those in low-risk group (<cut-off) ( Figure 1A). The similar results were found for OS, DSS, DFS, and PFS in both patients with colon and rectal cancer ( Figures 1B, C), indicating that prognostic significance of S100A4 mRNA expression does not depend on cancer type.
To investigate further, whether S100A4 mRNA expression could serve as an independent prognostic criterion of death or recurrence/progression, uni-and multivariate Cox regression analyses were applied. Univariate Cox regression analysis showed that S100A4 mRNA expression more than cut-off has prognostic value for poorer OS (HR=2,19; 95% CI   Table S6). In rectal cancer patients, COX analysis did not show the prognostic significance of S100A4 mRNA expression.
Statistical analysis showed that high S100A4 expression was associated with positive vascular invasion in RC patients (11,23 ± 1,26 vs. 10,75 ± 1,07; p=0,0163) ( Figure 1D). No significant associations were found with other clinical-pathological parameters. In CRC patients and colon cancer patients, no significant differences in S100A4 gene expression that was related to clinical-pathological parameters were found.  3.2 Elevated SPARC mRNA expression is an accurate independent prognostic factor for poor DFS and PFS in colon but not rectal cancer SPARC (secreted protein acidic and rich in cysteine, also known as osteonectin or BM-40) is a calcium-binding matricellular protein.
In the TME, SPARC is anti-angiogenic and affects tumor growth, extracellular matrix deposition (19). Similar to S100A4, ROC analysis allowed to divide all patients into high-and low-risk groups for prognosis of overall death, death from the disease, recurrence and progression based on SPARC mRNA level (Supplementary Table S7). SPARC mRNA expression more than cut-off in the high-risk group was significantly associated with worse DSS (p=0,0017) and worse PFS (p=0,0039), but with better OS (p=0,0401) in patients with CRC (Figure 2A). Similar tendency was shown for colon cancer patients ( Figure 2B). In rectal cancer, SPARC mRNA expression more than cut-off was associated with poor OS (p=0,0061), poor DFS (p=0,0371) and poor PFS (p=0,0120) ( Figure 2C) Table S10). Thus, we concluded that SPARC mRNA expression could serve as an independent prognostic factor for DFS and PFS in colon but not rectal cancer.

Increased SPP1 mRNA expression is an independent unfavorable criterion for PFS in both rectal and colon cancers
Secreted phosphoprotein 1 (SPP1, OPN, osteopontin) is an integrin-binding matricellular protein that has been found to be involved in many cellular processes such as cell signaling pathways, cell adhesion and migration, cell-mediated immunity, angiogenesis, and metastasis (20).
Univariate Cox regression analysis showed that mRNA levels of SPP1 more than cut-off were associated with decreased rates of DFS  (Supplementary Tables S12, S13).
Using CIBERSORT method, we demonstrated that mRNA expression of SPP1 and SPARC is significantly associated with M0, M1 and M2 macrophage phenotypes, naive B cells and CD4 T cells, while S100A4 gene expression correlates with M2 macrophage phenotype ( Figure 4B). SPP1 and SPARC mRNA expression strongly correlated with the expression of MSR1, CD163, MRC1 and MARCOmarkers associated with M2 TAM phenotype ( Figure 4B).
Additional Nanostring analysis allowed us to reveal that SPP1 is differentially expressed in immune CD45+ and tumor cytokeratin (CK)+ regions. The gene expression of SPP1 was higher in CD45+ compartments compared to CK+ regions (FDR=0.07, L2FC=0.66). SPP1 expression in the distinct regions of CRC was mostly associated with macrophages ( Figure 4C).
To confirm the expression of S100A4, SPP1, and SPARC in TAMs, we performed three-color IF analysis of human CRC tissue using confocal microscopy. IF analysis demonstrated that S100A4, SPP1 and SPARC were all expressed in CD68+ TAMs ( Figure 4D).
Thus, we demonstrated that levels of S100A4, SPP1 and SPARC are mostly inherent to the stromal component, in particular macrophages. Taken into account the observations made above, we used stromal-derived protein expression of S100A4, SPP1, and SPARC for further survival and correlation analysis.
3.5 Stromal levels of S100A4, SPARC and SPP1 retained an unfavorable prognostic value for patient outcome but became favorable for the pathological tumor parameters Protein levels of S100A4, SPARC and SPP1 remained unfavorable parameters for survival. In univariate COX analysis and Kaplan-Meier survival analysis, the protein level of S100A4 in high-risk group (expression more than cut-off meaning) remained prognostic only for OS in both combined CRC group of patients (HR=2,141; 95% CI [1,978], p=0,016) and patients with colon cancer (HR=2,679; 95% CI [1,136-6,320], p=0,024) ( Figures 5A, B). In multivariate Cox analysis, S100A4 mRNA expression was not an independent parameter. High protein expression of SPARC and SPP1 was associated with shorter recurrence-free survival (RFS) and PFS, respectively, in patients with rectal cancer (Figures 5C, D,  respectively). However, these parameters were not prognostic according to univariate and multivariate Cox analysis.
3.6 Neoadjuvant chemotherapy/ chemoradiotherapy reverse the prognostic value of S100A4 Finally, we found that neoadjuvant chemotherapy (NACT)/ chemoradiotherapy (NCRT) can reverse the activity of S100A4 from the pro-tumor to favorable one. In general, it can be hypothesized, that pro-angiogenic factors can induce formation of blood vessels with different functionalitywith different permeability for soluble factors or infiltration of immune cells from one side, and different permeability for cancer cells enhancing metastasis from another side (21,22). It can be assumed that specific type of vasculature can be beneficial for tumor growth before the treatment onset, while the same type of vasculature can be converted to the favorable for patients once chemotherapy is applied (23,24). The possible explanation of such effect can be enhanced permeability for the chemotherapy agent or for the selective anti-tumor immune cells. The idea of this study was based on our recent observation that chemotherapy can induce reprogramming of TAMs and launch the dysregulation in the expression of angiogenesis-associated factors (data not shown). Here, we found that stromal levels of S100A4 after treatment were lower in patients treated with NACT/NCRT compared to untreated patients (29,24 ± 22,49 vs. 19,04 ± 17,97, p=0,002), indicating that chemotherapy-based treatment can suppress its expression ( Figure 6A). Surprisingly, post-treatment stromal expression of S100A4 was higher in patients who have better response to NACT/NCRT (20,17 (12,87,35,00) for TRG1-2 vs. 12,78 (6,31;26,52) for TRG3-5, p=0,0438) ( Figure 6A). Response to NACT/NCRT was estimated by Mandard Tumor Regression Grading (TRG) system.
Our previous observations indicated that S100A4 can be expressed by both M1 and M2 macrophages in tumor tissues (unpublished data). We suppose that chemotherapy can induce repopulation of S100A4-expressed M1 and M2 macrophages in tumors.

Discussion
Currently, CRC prognosis is largely based on clinical and pathological parameters and focuses on the cancer stage at the time of diagnosis (25). Clinically validated prognostic biomarkers that can identify "high-risk" CRC patients are currently missing (25). Among cells of immune infiltrate, only T lymphocytes were included in Immunoscore classifier that was proposed for survival prediction: patients with "hot" tumors (where CD3+ and CD8+ T cells were detected) exhibit better RFS than patients with "cold" tumors (26,27).
Here, for the first time, we performed the complex analysis of mRNA expression using TCGA and GEO datasets, and protein expression using quantitative IHC of clinical samples for crucial regulators of tumor angiogenesis and tumor progression, expressed by tumor-associated macrophages: S100A4, SPP1 and SPARC (13). We considered colon and rectal cancer as two tumor entities as B C D A FIGURE 5 The association of stromal expression of S100A4, SPARC and SPP1 with survival and clinical-pathological parameters. (A) S100A4 stromal expression predicts poor OS for CRC patients and colon cancer patients, but is negatively correlated with tumor stage in rectal cancer. (B) The prognostic significance of S100A4 protein levels for overall survival in patients with CRC and CC revealed by univariate and multivariate COX analysis. (C) SPARC stromal expression is an unfavorable parameter for RFS in rectal cancer patients. SPARC protein expression is negatively associated with clinicalpathological parameters in combined CRC group and in rectal cancer patients. (D) SPP1 stromal expression is unfavorable for PFS in patients with rectal cancer. SPP1 protein expression is negatively associated with clinical-pathological parameters in combined CRC group and in rectal cancer patients. Log-rank test p-values are shown in Kaplan-Meier plots. Mann-Whitney U test was applied for the comparison of two groups. Box plots depict protein expression (min, Q1, median, Q3, max).
accumulating clinical data showed their differences in the prognosis and treatment strategies (28).
Using TCGA, survival analysis demonstrated that S100A4, SPP1 and SPARC could be promising predictors for the unfavorable outcome. In particular, S100A4 accurately predicted poor survival with high sensitivity and specificity for patients with CRC independently on cancer type, so S100A4 can be a universal unfavorable predictor in colon and rectal cancers. SPARC mRNA expression was a more specific marker for colon cancer patients. SPP1 mRNA expression was found to have a more significant predictive value for PFS in both rectal and colon cancers. Multivariate COX analysis also revealed that prognostic model consisting of S100A4 (HR=8,43; 95% CI [5,426], p<0,0001), SPARC (HR=1,86; 95% CI [1,141-3,041], p=0,012) and SPP1 (HR=1,86; 95% CI [1,124-308], p=0,058) can be established for the PFS in CRC cohort. Literature data indicate that in several cohorts of CRC, high SPP1 gene mRNA expression was associated with shorter OS, higher S100A4 expression with shorter DFS and OS, and high SPARС expressionwith worse DFS (29)(30)(31)(32)(33). These cohorts include patients with I-IV stages of CRC or only colon cancer, or only rectal cancer. We for the first time performed complex analysis of S100A4, SPP1 and SPARC mRNA levels in terms of survival rates in common CRC groups and in patients with colon and rectal cancer separately.
Next, we found significant correlation of S100A4, SPP1 and SPARC with the amount and M2 phenotype of TAMs. Confocal analysis confirmed the co-expression of these proteins in CD68+ TAMs in human CRC tissue. The co-expression of S100A4 and SPARC in CD68+ TAMs in human CRC tissue was demonstrated by us for the first time. We were able to find only one study describing the co-localization of CD68 and SPP1 in tumor stromal components in human CRC (34). Additionally, Nanostring technology allowed to find elevated expression of SPP1 in CD45+ immune compartment compared to CK+ tumor cell compartment as well as strong correlation of SPP1 with macrophage count in CRC tissues.
We found that, in contrast to mRNA expression level, S100A4, SPP1 and SPARC, stromal level negatively correlated with clinicalpathological parameters. These data correspond to the results found for other cohorts analyzed by IHC using tissue microarrays. In a cohort of 134 patients with CRC, negative correlations were found between SPP1 expression and distant metastasis, tumor invasion, tumor grade, and recurrence (35). In a cohort of 114 patients with colon cancer, a significant negative association was observed for SPARC expression in mesenchymal and stromal cells with the differentiation of tumors (36). For S100A4, several studies did not show any association of its protein expression with pathological parameters other than survival rates (32,37,38).
Interestingly, that stromal expression of SPARC and SPP1 was prognostic for the reduced RFS and PFS, respectively, only in rectal cancer patients, while S100A4 correlated with poor OS rates in CRC and colon cancer patients. Such controversial data can be explained by the dual role of TAMs in CRC progression. A few reports indicate that high amounts of TAMs that are the most abundant innate immune cell population in CRC are beneficial to CRC patients (39, 40).
Finally, we found that neoadjuvant chemotherapy (NACT)/ chemoradiotherapy (NCRT) can reverse the association of S100A4 with prognosis from the tumor progression-associated to favorable one. Accumulating evidence showed that chemotherapy can induce repolarization of macrophages in the TME (8). Chemotherapy-based treatment suppressed stromal expression of S100A4. Notably, S100A4 stromal levels were higher in patients with better response to NACT/ CRT, and S100A4 mRNA levels predicted better DFS among nonresponders. It can be explained by the mechanism of re-population of TAMs after chemotherapeutic intervention. It can be also assumed that B A FIGURE 6 S100A4 is associated with improved outcomes in patients undergone neoadjuvant chemotherapy-based treatment. (A) Stromal expression of S100A4 decreased in treated patients and correlated with favorable response to NACT/CRT. (B) S100A4 and SPARC mRNA expression decreased in nonresponders with progression. S100A4 high-risk score predicts better DFS. Mann-Whitney U test was applied. Box plots depict protein (A) and gene (B) expression (min, Q1, median, Q3, max). Log-rank test p-value is shown in Kaplan-Meier plot. specific type of vasculature can be beneficial for tumor growth before the treatment onset, while the same type of vasculature can be converted to the favorable for patients once chemotherapy is applied (23,24).
Thus, we demonstrated high prognostic significance of angiogenesis-associated factors S100A4, SPP1 and SPARC that are produced by TAMs in colorectal cancer. Our findings can help improve the prognosis of patients with CRC based on S100A4, SPP1 and SPARC expression levels. S100A4, SPP1 and SPARC can be helpful targets for developing novel immunotherapy and antiangiogenic therapy approaches.

Data availability statement
All data needed to evaluate the conclusions in the paper are presented in the paper and the Supplementary Material. Additional data related to this paper may be requested from authors. GEO data set was accessed via NCBI GEO Repositorium: https://www.ncbi.nlm. nih.gov/geo/query/acc.cgi?acc=GSE190826. TCGA-COAD and TCGA-READ data sets were obtained from NIH GDC data portal: https://portal.gdc.cancer.gov/. Results of NGS-GeoMx DSP analysis are presented in NCBI data portal (GSE221924): https://www.ncbi. nlm.nih.gov/geo/query/acc.cgi?acc=GSE221924.

Ethics statement
This study was performed in line with the principles of the Declaration of Helsinki. Approval was granted by the local committee of Medical Ethics of Tomsk Cancer Research Institute. Informed consent was obtained from all individual participants included in the study.

Author contributions
The study was conceptualized by IL and JK; EK, MR, TS, LT performed experiments and analyzed data; PI prepared the bioinformatics data analysis; AT, SA, AD, LZ and NC provided clinical data and enrolled the patients; IL, EK and MR wrote the manuscript with input from all authors. IL and JK interpreted data and contributed to the discussion.